function z=logdomain_diff(x,y)
%     x > y
%     z = lnc, y = lnb, x = lna
%     c = a - b
z = x + log(1 - exp(y-x));
end